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Accurate calculations of the fine and hyperflne splitting of the 2P state in Li and Be + isotopes using the 
explicitly correlated Hylleraas basis set are presented. Theoretical predictions including the mixing of P1/2 
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experimental values. It is concluded that precise spectroscopic determination of the nuclear magnetic moments 
requires elimination of nuclear structure effects by combining measurements for two different states. 
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I. INTRODUCTION 



The calculation of relativistic effects in the atomic struc- 
ture is most often performed with the explicit use of the Dirac 
equation, as in relativistic configuration interaction Qh, many 
body perturbation theory 0], relativistic coupled-cluster yj], 
or multi-configuration Dirac-Fock [4] methods. For light 
atomic systems the more accurate approach is based on the 
expansion of the energy in the fine structure constant a. This 
method allows for a systematical inclusion of relativistic and 
quantum electrodynamics (QED) contributions, as each cor- 
rection can be expressed in terms of the expectation value of 
some operator with the nonrelativistic wave function. With 
the use of explicitly correlated basis functions, the nonrela- 
tivistic Schrodinger equations for few electron systems can be 
solved very accurately. The high precision is achieved also for 
relativistic and QED corrections, provided more complicated 
integrals with inverse powers of inter-electronic distances can 
be performed. Such calculations, which rely on expansion in 
a, have been performed for hydrogen and hydrogen-like ions 
up to the very high order of ma 8 [5]. Slightly lower preci- 
sion was achieved for the helium fine structure and for other 
helium energy levels, all terms up to m a 6 order have been 
obtained with approximate inclusion of dominant m a 7 cor- 
rections 1 6]. For 3- and 4-electron atoms calculations have 
reached the order to a 5 with partial inclusion of to a 6 terms, 
which come from the electron self-energy. The complete cal- 
culation of the m a 6 contribution for 3-electron systems has 
not been performed so far. 

In this work we present accurate calculation of the fine 
and hyperflne splittings in Li and Be + ions through ma 4 
and to a 5 orders including the finite nuclear mass corrections. 
Lithium fine structure have already been calculated in Hyller- 
aas functions by Yan and Drake in |7[], but in a relatively small 
basis and with the neglect of P 1 / 2 and P 3 / 2 mixing, which we 
find to play a significant role in the isotope shift. The hyper- 
fine splitting of P-states was calculated in many works using 
explicitly relativistic methods yj, |2j, Li |4J] and with the non- 
relativistic multi-configuration Hartree-Fock method in HIH]. 
We find by a comparison with our results, that the most accu- 
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rate previous calculation was that performed by Yerokhin in 
QD. For the comparison with experimental values we include 
0(a 2 ) relativistic corrections from [1], known 0(a 2 ) QED 
corrections and draw a conclusion that a largest uncertainty 
comes from the not well known nuclear structure effects. 



II. FINE AND HYPERFINE OPERATORS 



Let us briefly start with the description of the fine and hy- 
perflne splitting in an arbitrary few electron atom. The fine 
structure, neglecting relativistic 0(a 2 ) corrections, can be ex- 
pressed as the expectation value with the nonrelativistic wave 
function of the following operator 



(.9-1) 



-5 — r a x p a 

1 £ 



to ton 



■ r a x Pn 



So [g r a b xpi- (.g - 1) f ab x p a ] 



(1) 



where g is the free electron g-factor, which includes here all 
QED corrections, Z is the nuclear charge in units of the ele- 
mentary charge e, to, ton are the electron and nuclear masses 
respectively, finally s a is the electron spin operator. For con- 
venience of further calculations we express Hf s in terms of F* 
and 4 elementary operators /* in atomic units, namely 
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where e — m a 4 and 



where h operators (in atomic units) are 
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The hyperfine structure, neglecting relativistic 0(a 2 ) cor- 
rections, is given by H^s operator in Eq. ©. We will treat 
the nucleus as any other particle with mass and with the 
g-factor g^ which is related to the magnetic moment /1 by the 
formula 
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where /in is the nuclear magneton and I is the nuclear 
spin. Nuclear masses, spins, magnetic dipole and electric 
quadrupole moments of Li and Be isotopes are taken from 
literature and are all presented in TableU With the help of <?n 
i?hfs can be written as 
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where Q is the electric quadrupole moment. For convenience 
of further calculations we express in terms of H a , , 
if*, and 77 y operators, namely 
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III. MATRIX ELEMENTS 

Matrix elements of the fine and hyperfine operators are 
evaluated with the nonrelativistic wave function. This func- 
tion is obtained by solving the Schrodinger equation in the 
3-electron Hylleraas basis set. Finite nuclear mass corrections 
are included by reduced mass scaling and perturbative treat- 
ment of the mass polarization correction. All matrix elements 
are expressed in terms of Hylleraas integrals, which are ob- 
tained with the help of recursion relations ITol [Til [T2I [l3ll . 
The high accuracy is achieved by the use of a large number of 
about 15000 Hylleraas functions, and we have already demon- 
strated the advantages of this approach by the calculation of 
the isotope shift in Li IU4I1 and Be + ions II 1511 . 

The nonrelativistic wave function is the antisymmetrized 
product of spacial and spin functions of the form 
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(21) 

X = [a(l)/3(2)-/3(l)a(2)]a(3), (22) 

where a z a(.) = a(.) and a z (3(.) = — /?(.). Matrix elements 
of each operator, after eliminating spin variables can take the 
standard form 
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and, what we call, the Fermi form 
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with the assumption that the norm is X^=i(*l*)s = L The 
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TABLE I: Data for Lithium and Beryllium isotopes. Atomic binding energy of — —7.281 au, Ebc = —14.669 au. The value for the 
quadrupole moment of 7 Be is a theoretical estimate flql . 
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matrix element of the fine structure Hamiltonian becomes 

(H ia )j = 



ijk 



and the fine splitting is 
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The matrix elements of the hyperfine structure Hamiltonian 
takes the form 
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where Aj and Bj are magnetic dipole and electric quadrupole 
hyperfine constants. They are all expressed in terms of stan- 
dard and Fermi matrix elements, namely 

Aj = T ,r.^ (J-G)j, (28) 
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Numerical values for all matrix elements involved in these cal- 
culations are presented in Table ITT1 They have been obtained 
by extrapolation to infinite basis set and uncertainties reflect 
the numerical convergence. Matrix elements of the fine struc- 
ture operators have been derived previously by Yan and Drake 
in fm and later by us in IU5I1 . Small differences with results 
of 0] come from the not very large number of basis functions 
used in that work. The hyperfine operators have been previ- 
ously obtained in several works, i.e. Illl9ll39[] and we compare 
our result with the most accurate one from [ 1 ] with which we 
agree well. 



IV. THE SECOND ORDER CONTRIBUTION 

The hyperfine Hamiltonian H\^ s mixes 2 2 Pi/ 2 with 2 2 P 3 / 2 
what leads to additional contributions to fine and hyperfine 
splittings 14011 . Since this mixing is not very large one can 
use the second order perturbative formula which involves off- 
diagonal matrix elements 
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To calculate them one can use Clebsch-Jordan coefficients and 
Racah algebra 04111 . In the simpler approach presented here, 
we introduce the operator K, such that (J, m\K\ J, to') = 
for J = 1/2, 3/2, but does not change L nor S, namely 
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Then the off-diagonal matrix elements can be transformed to 
the form 
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TABLE II: Matrix elements in atomic units of operators involved in the fine and hyperfine splitting of P-states, infinite mass and the mass 
polarization correction with the coefficient —m/(m + mjq). (k\h a \k) f corresponds to a c from Ref. Jl|], (i\h % J \j)f to 10 a s< j, e !:,fc Ws 
to 2 a;, and {i\h 1:l \j) s to b q /2, Numerical uncertainties are due to extrapolation to the infinite basis set and reflect the numerical convergence. 
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The second order correction to energy due to i/hfs in Eqs. (f30b 
neglecting the small Y 2 term becomes 
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where we omitted the magnetic octupole coupling, the so 
called Cj coefficient. Resulting corrections to the fine and 
hyperfine splittings are 
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V. RESULTS 



Numerical results for the fine splitting in Li and Be + iso- 
topes are shown in Table [TTTJ E^ is the leading contribu- 



tion with the exact electron g-factor, but in the infinite nuclear 
mass limit, E^' is the finite nuclear mass correction, and SEf s 
is the P\/2 — -F3/2 mixing term. The higher order relativis- 
tic and QED corrections are not known, as they have not yet 
been evaluated. Finally AE& is the isotope shift with respect 
to 7 Li and 9 Be + . Our result for this isotope shift in the fine 
structure AEf s of Li differs significantly from the previous 
calculations in J3l due to the inclusion of the important sec- 
ond order contribution 6Ef s . However, it differs also from all 
the experimental values, see Table Hill 

Numerical values of all significant contributions to the hy- 
perfine constants of the 2 2 Pf/ 2 and 2 2 P 3 / 2 states in Li and 
Be + isotopes are shown in Tables [IV] and [V] Ay 2 accord- 
ing to Eq. ([8]) involves the exact electron g-factor, and thus 
includes the leading QED corrections. The relativistic correc- 
tions A rel and B rcl have been calculated by Yerokhin in |l|] in 
terms of Gmi and Gei functions. Gmi is defined by 
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where relativistic corrections to Gmi are equal to 0.000015 
for 2 2 P 1/2 , -0.000039 for 2 2 P 3/2 states of Li, and 0.000153 
for 2 2 P 1/2 , -0.000161 for 2 2 P 3/2 states of Be+. These num- 
ber include also the so called negative-energy contributions. 
Ge2 is related to Bj coefficient by 
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where relativistic corrections to Gei for 2 2 P 3 / 2 are equal to 
-0.000 004 in Li and -0.000 013 in Be+. These relativis- 
tic corrections can in principle be evaluated within NRQED 
approach 15 111 but so far we have not been able to obtain ana- 
lytic formula for all Hylleraas integrals involved in matrix el- 
ements. The next to leading radiative (QED) correction A^ 
(beyond the anomalous magnetic moment) is proportional to 
the Fermi contact interaction and is known from hydrogenic 
atoms. In terms of the H a operator it is 
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TABLE III: Fine splitting of 2P-states in Li and Be+ isotopes in MHz with e = 2 Rca 2 = 6 579 683 921 MHz. AE ts is the isotope shift 
with respect to 7 Li and 9 Be. It is not clear whether the experimental value of Orth et al. 1421 for the 7 Li fine structure includes 8E{ S due to 
their diagonal and of-diagonal parametrization of hyperfine matrix elements. 
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TABLE IV: Hyperfine splitting of the 2P-states in Li isotopes in MHz. Results of Yerokhin Q|] are corrected by inclusion of 5 A and SB, and 
by the use of more accurate electric quadrupole moments for 6 Li and 7 Li. Results of Orth et al. I42H47I1 for A and B constants in 7 Li are 
shifted by 8 A and SB, as these authors parametrized results of their measurement by diagonal and of-diagonal parts separately. Uncertainties 
of final theoretical predictions are due to higher order corrections and the approximate treatment of the nuclear structure contribution. Not 
shown are uncertainties due to inaccuracies of magnetic dipole and electric quadrupole moments. 
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The last significant contribution is the finite nuclear size cor- 
rection, the extended electric and magnetic distribution within 
nucleus. It is given by the formula 
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rz = J d 3 rd 3 r'pE(r)p M (r') \r-r'\. (45) 
Using exponential parametrization of electric and magnetic 



formfactors 



Pe{t) 
Pm{t) 



3 2V 3, 
7rr| 



/r E 



3 V3 



-2VSr/r M 



(46) 
(47) 



1 Af 
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TABLE V: Hyperfine splitting of 2P-states in Be isotopes in MHz. 
Results of Yerokhin [l|] are corrected by inclusion of S A and SB. 
Uncertainties of final theoretical predictions are due to higher or- 
der corrections and the approximate treatment of the nuclear struc- 
ture contribution. Not shown are uncertainties due to inaccuracies of 
magnetic dipole and electric quadrupole moments. 

7 Be+ 9 Be+ "Be 4 " Ref. 

A T/2 -140.0696(3) -117.859 2(3) -504.8745(8) 

SA 1/2 -0.00961 -0.00683 -0.10520 

A *i/2 -0.096 8 -0.081 5 -0.349 MJ 



HI 

133] 
1501 

Ml 



Ml 

Ml 
Ml 



the Zemach radius rz is 

_ 35 (r E + r M ) 4 + 14 (r| - r 2 M f - (r E - r M ) 4 
rZ 32V3(r E +r M ) 3 

(48) 

For all but 11 Be nuclei we assume r E = thus 

r z = ^% = 1.263 r E , (49) 
16 v3 

and take charge radii from the recent isotope shift measure- 
ments in Li ll35ll an d Be + ll33ll supplemented with isotope shift 
calculations in IU5I1 . For the Gaussian distribution one obtains 
yj rz = 1.30 rg what demonstrates a weak dependence of 
rz on an arbitrarily assumed shape of the charge distribution, 
with one exception. The n Be nucleus has a single neutron 
halo, what means that r^j is much larger than r E and the nu- 
clear finite size becomes much larger. We employ here the 
result of direct calculations from 05211 . which is 

H l a as = H a (-0.000 717). (50) 

At the same time the nuclear polarizability correction is also 
much larger and of the opposite sign to the finite size effect. 
Since it is very difficult to estimate, it will be neglected 
here. The final results for ^1/2,^3/2, and B 3 / 2 include the 
uncertainty coming from the higher order corrections, which 
we estimate to be 25% of A qcd and the uncertainty due to 
the approximate treatment of the nuclear structure which we 



estimate to be 25% of A {ns for all A coefficients, while for 
the B coefficients we assume the final uncertainty to be the 
sum of 10%ofS rcl and 0.1% of (55. 

VI. CONCLUSIONS 

In comparison to experimental values we observe signifi- 
cant discrepancies for the isotope shift in the fine structure, 
see Table [Hi] Although the theoretical fine structure of 7 Li 
is consistent with experimental values, the differences can be 
associated to 0(a 2 ) relativistic corrections. The isotope shift, 
as it has already been noted in 10, HH] differs significantly be- 
tween different experiments and theoretical predictions. In 
view of the recent determination of the nuclear charge radii 
from the isotope shift of 2S , 1 / 2 — 2Pi/ 2 transition in Be + ions, 
it is important to resolve these discrepancies. In this respect, 
we note the recent critical examinations [53] of all experimen- 
tal values of the fine structure and isotope shift measurements 
in 6 Li and 7 Li. Considering hyperfine splittings we observe 
good agreement with the previous calculations of Yerokhin in 
G]], particularly for the A\i<i coefficients. Slight discrepancies 
with experiments for the A coefficients of the 2P state indi- 
cate that the magnetic moment obtained from the hfs measure- 
ment for the 25 state may not be as accurate as claimed. This 
is because the treatment of the nuclear structure corrections by 
the elastic charge and magnetic formfactors is very approxi- 
mate, and the accuracy of this approximation is not known. 
We think that the more accurate approach shall employ the 
effective nuclear Hamiltonian using the so called chiral per- 
turbation theory. Then the nuclear structure correction to the 
atomic hfs consists of the leading Low correction, Zemach 
corrections from individual nucleons and the nuclear vector 
polarizability 15411 . Unfortunately, the explicit calculations for 
nuclei with more than 3 nucleons is difficult and has not been 
performed so far. Certainly the nuclear vector polarizability 
correction is significant for halo nuclei, and it would be worth 
to calculate it. At present, without detailed knowledge of nu- 
clear structure, the determination of magnetic moments from 
atomic spectroscopy measurements can be uncertain. There- 
fore, better accuracy can be achieved when two measurements 
are combined in such a way, that this nuclear structure cor- 
rection, proportional to the Fermi interaction cancels out, for 
example in Ai/ 2 + A 3 / 2 of the P state of Li and Be + . Theoret- 
ical accuracy for this combination is limited only by higher or- 
der QED corrections and knowing both A constants, we shall 
be able to derive magnetic moments with relative precision of 
about 10~ 5 without referencing to magnetic moments of sta- 
ble isotopes, or with precision of the magnetic moment of the 
reference nucleus. 
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SA 1/2 

A l/2 
.qod 
^1/2 

4 tn , 3 

^1/2 
A l/2 



expt. 
expt. 



SA 



"■-3/2 

3/2 
ATel 
A 3/2 
.qod 
A 3/2 

4 fns 
^3/2 

A 3/2 



B 3.)°2 



SB 



3/2 
3 3/2 



3 3/2 



7 Be+ 
-140.069 6(3) 
-0.00961 
-0.096 8 
0.008 26 
0.010 85 
-140.157(3) 

-140.17(18) 

-1.21533(2) 
-0.003 90 

0.020 3 
-0.008 26 
-0.01085 
-1.218(3) 



-2.63619(1) 
-0.025 43 
0.000 20 
-2.66142(3) 



9 Be+ 
-117.8592(3) 
-0.006 83 
-0.081 5 
0.006 95 
0.008 69 
-117.932(3) 

-117.926(4) 
-118.00(4) 
118.6(36) 

-1.02481(2) 
-0.002 76 

0.0171 
-0.006 95 
-0.008 69 
-1.026(3) 
-1.018(3) 

-2.281 54(1) 
-0.018 03 
0.00017 
-2.299 40(3) 
-2.299 25(17) 



11 Be+ 
-504.874 5(8) 
-0.105 20 
-0.349 
0.029 75 
0.055 50 
-505.245(16) 

-505.41(5) 

-4.395 48(8) 
-0.052 60 

0.073 5 
-0.029 75 
-0.055 50 
-4.460(16) 
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